Universal algebraic relaxation of fronts propagating into an unstable state 
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We analyze "pulled" or "linearly marginally stable" fronts propagating into unstable states. 
While "pushed" fronts into meta- and unstable states relax exponentially, pulled fronts relax alge- 
braically, and simultaneously the standard derivation of effective interface equations breaks down. 
We calculate all universal relaxation terms of uniformly translating pulled fronts. The leading 1/t 
and l/t s ^ 2 corrections to the velocity are determined by the dispersion relation of the linearized 
equation only. Our analysis sheds new light on the propagation mechanism of pulled fronts. 

PACS numbers: 02.30.Jr, 47.54.+r, 47.20.Ky, 03.40. Kf. 
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Consider systems far from equilibrium with a contin- 
uous order parameter, where a stable state propagates 
into an unstable state, and assume, that thermal pertur- 
bations can be neglected. Some experimental examples 
we discuss below. If the initial profile is steep enough, 
arising, e.g., through a local initial perturbation, it is 
known that the propagating front in practice al- 

ways relaxes to a unique profile and velocity. Depending 
on the nonlinearities, one can distinguish two regimes: as 
a rule, fronts whose propagation is driven ("pushed") by 
the nonlinearities, resemble very much fronts propagating 
into metastable states and their relaxation is exponential 
in time. This regime is often refered to as "pushed" 
or "nonlinear marginal stability" jl. If, on the other 
hand, nonlinearities mainly cause saturation, fronts prop- 
agate with a velocity determined by linearization about 
the unstable state, as if they are "pulled" by the linear in- 
stability ( "pulled" or "linear marginal stability" regime) . 
In the simplest case [Eq. (fj])] of a pulled front it has 
been proven rigorously Q that the velocity relaxation is 
proportional to 1/t in leading order, and there are some 
heuristic arguments Q, that this is generally true for 
pulled fronts. In this paper, we identify the general mech- 
anism leading to slow relaxation of uniformly translating 
fronts and use it to introduce a systematic analysis which 
allows us to determine all universal asymptotic terms. 
We get a unification and extension of various seemingly 
unrelated approaches 0-0] as a bonus. 



Our present investigation actually was motivated by an 
intimately related problem, namely the derivation of ef- 
fective interface equations for some order parameter field 
varying on a short length scale and coupled to some ex- 
ternal field (like temperature in a solidification or com- 
bustion front) varying on some outer length scale |)|. 
In such problems, it is well known that in the limit in 
which the curvature of the front or domain wall is much 
smaller than its width, smooth interface problems reduce 
to those with a sharp interface with boundary conditions 
that are local in space and time. This connection actu- 
ally lies at the basis of recent advances in the numerical 
studies of interfacial growth problems like dendrites po) . 



Although this is often not made explicit, this reduction 
of a smooth front problem to a sharp interface formula- 
tion amounts to an adiabatic decoupling of the interface 
relaxation from the dynamics of the outer field (s). This 
decoupling is only possible if the profile of the front prop- 
agating into a metastable state relaxes exponentially fast 
to its asymptotic shape and speed. If the order parameter 
front is a pulled front, the standard derivation of effec- 
tive interfaces (^] breaks down and at the same time the 
relaxation becomes non-exponential. Both are due to the 
properties of the spectrum of the operator C* derived by 
linearizing about the asymptotic front, as we will discuss 
in more detail elsewhere jll]] . 

We here analyze a pulled uniformly translating front 
in one dimension without external fields and its algebraic 
relaxation towards its asymptotic shape. We find that 
not only every sufficiently steep initial profile relaxes to 
a unique asymptotic front profile, but that also this pro- 
file is asymptotically approached along a unique trajec- 
tory. Such properties, that are independent of the precise 
initial conditions, we call universal. We construct all uni- 
versal terms explicitly in an asymptotic 1 / \ft expansion 
and confirm our predictions numerically. We not only 
calculate the dominant 1 jt term of the velocity HQ , but 
also find the subdominant 1/t 3 / 2 term. The term 1/t 2 is 
nonuniversal, since it at least partially can be generated 
from the 1/t term through a temporal shift of the initial 
conditions t — > t + tQ. Furthermore we, for the first time, 
predict the relaxation of the shape of the profile. Techni- 
cally, we here construct the contracting center manifold 
in function space about the asymptotic front profile. 

We use the word "universality" in analogy to scaling 
in critical phenomena: Not only do a large class of ini- 
tial conditions reach a unique asymptotic profile (fixed 
point) in an infinite-dimensional function space, but also 
the asymptotic approach towards this fixed point is com- 
pletely determined by the expansion about it. This is 
reminiscent of the universal corrections to scaling. 

Fronts propagating into an unstable state arise in var- 
ious fields of physics: they are important in many con- 
vective instabilities in fluid dynamics such as the onset 
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of von Karman vortex generation |L2l , in Taylor Ej|] and 
Rayleigh-Benard flTij convection, they play a role in spin- 
odal decomposition near a wall Eg], the pearling instabil- 
ity of laser- tweezed membranes Eg], the formation of ki- 
netic, transient microstructures in structural phase tran- 
sitions [jl7|, the propagation of a superconducting front 
into an unstable normal metal Eg], or in error propaga- 
tion in extended chaotic systems ]19[ . Our own interest in 
the relaxation of pulled fronts stems from our search for 
an interfacial description for so-called streamers, which 
are dielectric breakdown fronts poj . The experimental 
relevance of front relaxation is illustrated on propagating 
Taylor vortex fronts. Here the measured velocities were 
about 40% lower than predicted theoretically and only 
later numerical simulations |2l|] showed that this was due 
to slow transients. 

Our analysis can be formulated quite generally for par- 
tial differential equations which are of first order in time 
but of arbitrary order in space, as long as they admit 
uniformly translating "pulled" solutions. ("Pulling" is 
explained in the next paragraph.) For ease of presenta- 
tion we guide our discussion along two examples which 
we have investigated analytically as well as numerically. 



d t ct>(x,t)=d*<l> + f(<P) , /(</,) =0-0 3 
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with cj>, x, t real, is the nonlinear diffusion equation, ft is 
a prototype equation for pulled fronts. It is also known 
as KPP equation (after Kolmogorov et ah), Fisher equa- 
tion, or FK equation. In (]]]), the state 4> = is unstable 
and the states <fi = ±1 are stable. We consider a situation 
where initially 4>(x, 0) asymptotically decays quicker than 
e~ x for large x, or in particular one with <f>(x, 0) ^ in 
a localized region only. The region with <fi ^ expands 
in time, and a propagating front evolves. It has been 
proven rigorously, that relaxation is always to a unique 
front profile (f>*(x — v*t) with velocity v* — 2 0, and that 
the velocity relaxes asymptotically as v(t) = 2 — 3/(2i) 
H . The extension 



d t 4>(x, t) = d 2 x 4> - 7^0 + <t> - 4? 



(2) 



is often refered to as the EFK (extended FK) equation, 
and serves as a model equation for higher order equations 
that admit front solutions. The rigorous methods of [QJ|] 
are not applicable here, but the algebraic relaxation be- 
havior towards a unique uniformly translating front with 
velocity v* continues for < 7 < 1/12 |2^fl. 

Since the basic state <fi = into which the front prop- 
agates, is linearly unstable, even a small perturbation 
around </> = grows and spreads by itself. According 
to the linearized equations any localized small perturba- 
tion will spread asymptotically for large times with the 
"linear marginal stability" speed v* Q] 



d^SLU 



<93 k 



- v* = 



d^k 



= 



k> 



%Uj{k* 



= v* , (3) 



where Lo(k) is the dispersion relation of a Fourier mode 
e -iut+%kx _ rj,^ £ rg ^. ^ wQ e q ua tj ons j n (||) are saddle point 

equations in the complex k plane that govern the long 
time asymptotics of the Green's function in a frame mov- 
ing with the leading edge of the front. The third equa- 
tion expresses that for selfconsistency, the linear part of 
the front should neither grow nor decay in the comoving 
frame. If the nonlinearity pushes the front to a velocity 
ip > v* , we call the front nonlinearly marginally stable or 
pushed. If the front propagates with velocity v* , we call 
the front linearly marginally stable or pulled. Replacing 
the nonlinearity in (|l]) by e.g., /(</>) = ecj) + (j) 3 — </> 5 , 
one finds the fronts to be pushed for e < 3/4, and to be 
pulled for e > 3/4 |J]. We stress that for a given partial 
differential equation, the pulled speed v* and the asymp- 
totic shape of the profile are determined explicitly by the 
dispersion relation ||^] . 

Our analysis and predictions can now be summarized 
as follows: for Eqs. (§) and (§) ( with < 7 < 1/12), 
and others where localized initial conditions develop into 
uniformly traveling pulled fronts, i.e., into fronts with 
asymptotic speed v* and = = 5Rcj(fc*), the asymp- 
totic relaxation for t — * 00 is given by 
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v*t-X(t), (4) 
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Here D = ±d 2 %uj/ (d%k) 2 \ k * plays the role of a diffu- 
sion coefficient for perturbations about the asymptotic 
shape, and A = > is the asymptotic decay rate of 
(j)*: (f>*(£) oc (£ + const.) e~ A? for £ -> 00. The factor 
(£ + const.) is due to two roots of w{k) — vk coinciding, 
which is implied by the saddle point equations (||) . Note 
that all quantities in (|^) are determined by the linear 
dispersion relation oj(k). 

It is central to the analysis, which we will further elu- 
cidate below, that the coordinate system, in which the 
asymptotic shape 4>* is subtracted, moves with the speed 
v* + X of 0(£, t). Only then the shape correction n stays 
small for all times. We then find through some expansion 
in the "interior region" of the front, where \r]\ <g; 0*, 

€^t + x(0 + or 2 ) (6) 

= + Xn sh (0 + 0(t- 2 ) , 7 hh = (5(f> v /5v) . 

V* 

In the far edge, where £ > 0(V Dt) 1, a different ex- 
pansion is needed, as the transient profile <j) falls off faster 
than 0*, so that 77 « — (f>*. Matching to the interior (|^) 
and imposing that the asymptotic shape (f>* is approached 
for t — ► 00, and that the transients are steeper than e _A ^ 
for £ — > 00, uniquely determines the velocity correction 
X (g) and the intermediate asymptotics 

</>(£, t) » e - A «-« 2 /(4^) (£ + cons t. + 0(1/ Vi)) . (7) 
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Let us now first discuss three physical observations, 
that are crucial ingredients of our systematic analysis. 
We then present our numerical results and finally briefly 
sketch our calculation. 

(i) The leading 1/t term in X(t) in (|) can be under- 
stood intuitively through a heuristic argument [Q: If we 
start from localized initial conditions and analyze in the 
asymptotic frame £* — x — v*t, ,t) should approach 
<$>*({;*) as t — > oo, but for a fixed time, <j> should fall 
off faster than <p* as £* —> oo. To study this crossover, 
consider for simplicity Eq. ([!]); if we linearize, and sub- 
stitute 0(r , t) = e-^Vlf » t) (with v* = 2, A = 1 = D 
in this case), we get the simple diffusion equation dt4> = 
d 2 ,ip. Clearly, the similarity solution which matches to 

r(D~e- A «*(r+ const.) is V ~ (C A 3/2 )e"«* 2 / 4 *, so 
^el-Ar-s^int-Hnfr-r 2 /^] ||. Hence, if we now track 
the position £? of constant height /i<Cl, which is defined 
as </>(£*, t) = h, we find = -3/(2A) hit + . . . in the 
frame £*. This is precisely the leading term of -X"(t)! 

(m) Observation (i) shows that within the frame £* = 
x — v*t, the leading edge of the profile moves back a dis- 
tance — 3/(2A)lni. But this must be true for all heights 
h, not just for h <C 1: The front width W is finite in 
equations like (Q) and (g). So positions are re- 

lated through - = 0(W) = 0(1). In other 

words, the dominant shift — 3/(2A)lnt determined from 
the leading edge h <C 1, is the same for all amplitudes 
h. This is why in order to study the asymptotics, we 
have to perturb about the asymptotic shape 0*(£) and 
not about </>*(£*), since the two profiles in frames £ or 
£* are pulled arbitrarily far apart due to the diverging 
logarithmic shift X(t). A very unusual situation indeed! 

(Hi) In (g), 77 s /j is the shape mode 8<f> v /5v\ v * which 
gives, to linear order, the change in the shape of the uni- 
formly translating profile <p v when varying v about v* 
p4j . Since the instantaneous velocity of the profile po- 
sition is just v* + X(t) (< v*), (||) expresses that up 
to order t~ 2 the front profile in the interior front region 
is just the uniformly translating profile 4> v with velocity 
v — v* + X(t). This was actually conjectured for Eq. (|l|) 
by Powell et al. g , based on numerical observations, and 
our analysis gives its first derivation in leading order; in 
order 1/t 2 the conjecture does not hold. 

We have tested our predictions by numerically inte- 
grating Eqs. (Q) and (|^) forward in time, starting from lo- 
calized initial conditions. In Fig. 1(a) and (b), we present 
velocities Vh{t) of various amplitudes h, in (a) for Eq. ([!]), 
and in (b) for Eq. (§]) with 7 = 0.08. Note that the crit- 
ical value of 7 is 7c = 1/12 = 0.083. Eqs. @-@ imply 
that in the lab frame x, v h (t) = v* + X(t) + g(h)/t 2 
(where g(h) can be expressed in terms of r) s h and d^<p*). 
Thus we plot Vh(t) — v* — X(t) versus 1/t 2 for various h. 
According to our prediction, all curves should then con- 
verge linearly to zero as 1/t 2 — > 0. Clearly, the numerical 
simulations fully confirm this for both equations. 



For our prediction (|J) of the shape relaxation, the most 
direct test is to plot [</>(£, t) - (j)*(()}/(X(t)r] sh (^)) as a 
function of £ for various times. This ratio should con- 
verge to 1 for large times. As Fig. 1(c) shows, this is 
fully borne out by our simulations of the nonlinear diffu- 
sion equation (Q). Moreover, the crossover for large pos- 
itive £ is fully in accord with our result that the proper 
similarity variable in the far edge is £ 2 /t — see Eq. (Q). 

We finally give a brief sketch of the systematic anal- 
ysis, taking the nonlinear diffusion equation (Q) as an 
example. Full details will be published elsewhere pT| . 

We first consider the "front interior" region, where the 
deviation r](£,t) of <j> about </>*(£) is small, i.e., \rj\ <C </>*. 
As there is some freedom in choosing £ due to transla- 
tion invariance, we choose quite arbitrarily the condition, 
that <f>(0,t) = h = 4>*(0), so that 77(0, t) = 0, as was also 
done in Fig. 1(c). Substituting (g) into ([[]), we obtain 



d t rj — C*rj + X d^(r]- 



^fl r, 2 + 0(rf) , (8) 



£*=a$+v*d 6 +/V(0) • 



(9) 



The inhomogeneity Xd^cp* in (|8|) is due to the fact 
that 4>*(0 is a solution of (|l|) only if X = 0. Since 
X(t) = 0(t _1 ), and since in the front interior \t)\ <C <fi* , 
the inhomogeneity induces an ordering in powers of 1/t, 
which suggests an asymptotic expansion as 



x = c ± 

t 



C3/2 

t 3/2 



£2 
t 2 



(t *\ ^ J V3 / 2 



(10) 
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The necessity for actually expanding in powers of l/yi 
emerges from matching to the similarity solutions in the 
far edge. Substitution of the above expansions in (||) 
yields a hierarchy of o.d.e.'s of second order 

C*r]! = -c^cjf , C*r] 3/2 = -c 3/2 9 5 0* (12) 
C*7] 2 = -c 2 d^(t>* - cidfqi — 771 — f"((j)*)rj 2 l /2 etc. 

The hierarchy is such that the equations can be solved 
order by order. Each rji is uniquely determined by its dif- 
ferential equation, the appropriate boundary conditions 
and the requirement r]i(0) = 0. The equations for rji/ci 
resp. 773/2/C3/2 are precisely the differential equation for 
i] sh = 6<f> v /6v\ v *, cf. Eqs. (£§), (||) and paragraph (Hi). 

By expanding the rji for large £, one finds that they all 
behave like e~ A ^ = times a polynomial in £, whose 
degree grows with i. The r\i expansion is therefore not 
properly ordered for large £. This just reflects the fact 
that on the far right, 77 and <f>* must almost cancel each 
other. This is required for fronts that emerge from lo- 
calized initial conditions, whose total profile thus decays 
faster than tfi*. A detailed investigation of this region 
shows that z = ^ 2 /it is a proper similarity variable here, 
and suggests that here the proper expansion is 
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Vtg^.(z)+g (z) + ^- + 



(13) 



Upon substitution of this expansion into the original par- 
tial differential equation, linearized about = 0, we now 
find a different hierarchy of ordinary differential equa- 
tions for the functions g n /2( z )- I n this case, the condi- 
tions to be imposed on the g n /2^> is that they do not 
diverge as e z as z — > oo, and that they match, in the 
language of matched asymptotic expansions, the large £ 
"outer" expansion of the "inner" solution based on the 
rji [25|. These conditions fix the parameters c\ and C3/2 
in (It), and this yields the solution given in Eqs. (H)-(@) 
Jil"| | . The structure of the analysis is essentially the same 
for higher order equations like (|2|). 

In summary, our results show that the 1 jt relaxation of 
pulled fronts is essentially due to the crossover to a Gaus- 
sian shaped tip in the leading edge of the front. The non- 
linearities dictate the asymptotic tip shape (p* cx £e~ A ? 
for t — > 00 and £ large. This asymptote determines the 
coefficients and the 1 /i 3 / 2 term in the velocity correction 
X (|J). We finally note that analytical arguments as well 
as numerical simulations indicate that many of the above 
arguments can be generalized to the case of pattern form- 
ing fronts, occuring, e.g., in Eq. (§ for 7 > 1/12 or m 
the Swift-Hohenberg equation Q|. Work is in progress. 

This work was started in collaboration with C. Caroli 
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FIG. 1. (a) and (b): Velocity correction Avh(t) = Vh(t) — 
v h = x h and for t > 20. (a): Eq. (1), thus 7 = 0, A = 1 = D, v* = 2. (b): Eq. (2) with 7 = 0.08, thus D = 0.2, A = 1.29, 
v* = 1.89. (c) Data from (a) plotted as t) — <f>* {£,))/ (Xri a h(£)) over £ f° r various t. <f)*{V) (dashed) for comparison. 
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